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ABSTRACT 

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. 
In total 20 lines are detected, including CO J = 4 — 3 through J = 13 — 12, 6 H2O rotational lines, 
and [C 1] and [Nil] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in 
NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV 
and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We 
applied both C and J shock models to the H2 v = 1 — 5(1) and v = 2 — 1 S(l) lines and the CO 
rotational ladder. The CO ladder is best reproduced by a model with shock velocity v s = 10 kms -1 
and a pre-shock density ijh = 5 x 10 4 cm~ 3 . We find that the solution best fitting the H2 lines 
is degenerate: The shock velocities and number densities range between v s = 17 — 47 kms -1 and 
7Jh = 10 7 — 5 x 10 4 cm -3 , respectively. The H2 lines thus need a much more powerful shock than 
the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 kms -1 ) 
shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. 
This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO 
line-to-continuum ratio is a key diagnostic for the presence of shocks. 

Subject headings: galaxies: individual (NGC 6240) — galaxies: active — galaxies: nuclei — galaxies: 
starburst — infrared: galaxies 



1. INTRODUCTION 

We present HerscheR SPIRE FTS (IGriffin et aLllMlOft 
observations of the nearby luminous infrared galaxy 
NGC 6240 (IRAS 16504+0228, UGC 10592). With a 
redshift of z = 0.024576 and a 7-year WMAP flat cos- 
mology (H = 70.50 kms^Mpc- 1 , fi mat ter = 0.27, 
^vacuum = 0.73) NGC 6240 is at a luminosity distance of 
Z?l = 107 Mpc, with 1" corresponding to 492 pc. The de- 

1 Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with 
important participation from NASA 



rived 8 — 1000 /im luminosity of the merger galaxy NGC 
6240 is Lir = 7.5 x 10 11 L©. Different power sources are 
suggested for this infrared luminosity, which complicate 
the interpretation of observations. The source is a strong 
X-ray emitt er. Based on modelin g of Chandra between 
0.5 — 8 keV. IKomossa et all (|2003[ ) derive extinction cor- 
rected luminosities of Lx,n = 1-9 x 10 42 erg s~ x and 
Lx,s — 0.7 x 10 42 erg s _1 for the northern and southern 
AGN cores in the energy rang e 0.1 — 10 keV. How ever, 
using BeppoSAX observations. I Vignati et all (|1999h find 
that the intrinsic hard X-ray luminosity shows up above 
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9 — 10 keV, and their models of the higher energy yields 
L(2 - 10 keV) » 3.6 x 10 44 erg s" 1 . Hydrogen re- 
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nuclei indicate recent star formati on. The H2 v = 
1 — S (l) 2.12 ^m line, presented by Ivan der Werf et al.l 
(1993), shows a peak in the overlap region between these 
two nuclei, which are 2" apart. This H2 emission extends 
over several kpc and shows a complex morphology. The 
authors conclude that the bright H2 emission between the 
nuclei is generated in shocks resulting from the collision 
of the interste llar media (ISM) o f the merging galaxies. 
More recently, lEngel et al.l (|2010D observed the same line 
at high spatial (0.5") and spectral (~ 90 km s _1 ) reso- 
lution. The line profiles also indicate multiple compo- 
nents, and the dispersion map suggests that the gas is 
highly disturbed and turbulent. Interferometric observa- 
tions of several CO, HCN, and HCO + lines show a high- 
density peak between t he two nuclei (ITacconi et al.llT999l; 
INakanishi et~aTl 120051; llono et al.l [2007D . iTacconi et al.l 
(1999) conclude that NGC 6240 is in an earlier merging 
stage than Arp 220, the prototypical ultra-luminous in- 
frared galaxy, and state that the gas is in the process of 
settling between the two nuclei, and di ssipating angular 
momentum rapidly. lEngel et al.l (|2010t ) also presented a 
12 COJ = 2 — 1 interferometric map, and conclude that 
the H2 and CO emission are coextensive, but do not co- 
incide with the stellar emission distribution. 

This paper is ordered as follows. In Section 2, the ob- 
servations, data reduction and line luminosities are dis- 
cussed. In Section 3, we discuss the different excitation 
components (PDR, XDR, shocks) in combination with 
the geometry of NGC 6240. Consequently, we compose 
a CO ladder using both the available SPIRE FTS and 
ground based CO data from IPapadopoulos et al.l ()2011[ ) 
and compare this system to the UL IRG Mrk 231 (pre- 
vious ly studied by us with Herschel, Ivan der Werf et all 
201 01). Then we analyze this with shock models from 
Kristensen et al.l ([2007). In Section 4, we conclude with 
a discussion and a summary of the results. 

2. OBSERVATIONS, DATA REDUCTION, AND RESULTS 

NGC 6240 was observed (Observations ID 1342214831) 
in staring mode with the SPIRE FTS on February 27, 
2011, as part of the Herschel Open Time Key Program 
HerCULES (P.I. Van der Werf). The high spectral res- 
olution mode was used, yielding a resolution of 1.2 GHz 
over both observing bands: the low frequency band cov- 
ering v = 447 - 989 GHz (A = 671 - 303 fim) and 
the high frequency band covering v = 958 — 1545 GHz 
(A = 313 - 194 ^m). In total 97 repetitions (194 FTS 
scans) were carried out, yielding an on source integration 
time of 12920 s (3.6 hrs) for each band. A reference mea- 
surement was used to subtract the combined emission 
from the sky, telescope and instrument. The data were 
processed and calibrated using HIPE vers ion 6.0. The 
exten t of the CO J = 3 — 2 emission is 4" (jWilson et alJ 
120081 ). while the SPIRE beam varies from 17" to 42" over 
our spectrum. Therefore, a point source calibration pro- 
cedure was adopted, and no corrections for wavelength 
dependent beam coupling factors were necessary. 

The full SPIRE FTS spectrum of NGC 6240 is shown 
in Fig. [TJ In the overlap region between the two fre- 



quency bands (958 — 989 GHz), the noisy parts of the 
two spectrometer bands were clipped and plotted on top 
of each other. A total of 20 lines were detected of which 
one line at the observed frequency ^ rcs t = 1481.6 GHz re- 
mains unidentified. We also see hints of absorption and 
emission between, e.g., 1300 and 1350 GHz. The signifi- 
cance of these features are still under investigation, and 
requires a more advanced flux extraction than currently 
adopted. The identified lines include the CO J = 4 — 3 to 
J = 13- 12, 6 H 2 lines, [Ci] 370 and 609 fim, and [Nil] 
205 /im. A well-constrained upper limit was obtained for 
the 13 CO J = 6 — 5 line. The line fluxes are listed in 
Table [1] In all cases a simple integration over the line 
was done. This procedure was repeated with approxi- 
mately 100 different choices of baseline and integration 
intervals. Specifically, the integration is done from ±5 to 
±9 GHz centered on the line with steps of 0.4 GHz (10 
integration borders). For each of these integration bor- 
ders, the baseline window is set from ±2 to ±3.5 GHz 
with steps of 0.15 Ghz (10 iterations for each of the in- 
tegration borders). The average is taken as the line flux 
and the standard deviation of the different realizations is 
taken as the error, and excludes the relative calibration 
error of 7 percent. This is also done for fluxes of lines 
that are blended. For these, for each iteration, a double 
Gaussian was fitted and the peaks were used to set the 
fractional flux for each of the lines. As additional model- 
ing constrai nts, the table lists fluxes f or three lower CO 
transitions (jPapadopoulos et al.ll2011[). an d luminosities 
of four H 2 tr ansitions flTecza et alJ l2000( ) . From these 
four H 2 lines, iTecza et all 72000) estimated an unatten- 
uated power of Lh 2 « 2 x 10 9 L q . 

3. ANALYSIS 

We will focus on the excitation of the CO ladder, and 
compare CO excitation conditions to those of H 2 . A 
multicomponent Large Velocity Gradient (LVG) analy- 
sis of the CO ladder, considering also constraints from 
the HCN and [CI] lines will be presented in another 
paper (Papadopoulos et al., in prep.). Preliminary 
mass est imates from this study are c onsistent with thos e 
found bvFapadopoulos et all (l2012D.lGreve et al.l (l2009D. 
ITacconi et al.l (|1999f ). and lEngel et alJ (|2010D . Estimates 
range between M tota i = 3 x 10 9 and 4 x 10 10 M , de- 
pending on whether low and/or high density tracers are 
used. 

The total luminosity in the 12 CO lines listed in Table 
[I] is Leo ~ 5 x 10 8 Lq. In the absence of J up > 13 mea- 
surements this is a lower limit to the total luminosity in 
the CO lines. Although the line luminosity of the CO 
J = 8 — 7 transition is largest, the luminosity in the indi- 
vidual CO lines is only slowly decreasing for higher rota- 
tional quantum numbers, and a significant contribution 
is expected from transitions with J up > 13. The shape of 
the CO ladder of NGC 6240 is similar to that of Mrk 231 
(see Fig. ^ and can be fitted by two photon-dominated 
region (PDR) mo dels and an X-ray domi nated region 
(XDR) model (see Ivan der Werf et al.ll2010l ). The physi- 
cal and geometrical properties of NGC 6240 are different 
from those of Mrk 231. We argue below why the fit used 
for Mrk 231 is not appropriate for NGC 6240, and why 
shocks must be responsible for the CO excitation in NGC 
6240: 

Absence of OH+ and H 2 0+ : The NGC 6240 spec- 
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Fig. 1. — Full SPIRE FTS spectrum of NGC 6240 and zoom-in on line blends, with on the x-axis the frequency and on the y-axis the flux 
density in Jy. The detected lines are marked: CO (red), [Cl] and [Nil] (green), H2O (orange), and unidentified (brown). In the overlap 
region between the two frequency bands (958 — 989 GHz), the noisy parts were clipped and plotted on top of each other. 
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Fig. 2.— Compariso n of the CO ladder of N GC 6240 to those 
obtained for Mrk 231 <van der Werf et al.l l2010). The CO ladder 
of Mrk 231 is normalized to the CO J = 8 - 7 luminosity of NGC 
6240. 

trum does not show emission line features of the 
ionic species OH + and H2 O + , observed in Mrk 231 
(jvan der Werf et al.ll20loh . Large OH+ and H 2 0+ abun- 
dances are only sustained in gas clouds with high ion- 
ization fractions (x e > 1CP 3 ), which a re produced by 
eleva ted cosmic ray or X-ray fluxes (cf., iMeiicrin k et al.l 
I2011D . Their absence hints that the bulk of the gas is not 
exposed to high ionization rates resulting from AGN or 
starburst /supervae activity. 

Line-to- continuum ratio: The CO luminosity (up to 
the J = 13—12 transition) to infrared (measured between 
8 — 1000 /im) luminosity ratio is Lco/^m = 7 x 10~ 4 
in NGC 6240. This is exceptionally high, and approx- 
imately an order of magnitude higher than the ratio 



of NGC 6240 (Ivan der Werf et al.l 19931) . Our PDR and 
XDR models ([Meiierink fc Spaansll2005l : iMciierin k et al.l 
2007) give a maximum Lco/iiR ratio of ~ 10 4 (as- 
suming /(FIR) = 2.6 x 10 _4 Go erg s" 1 cm" 2 sr -1 , see 
Kaufman et al.l fl999f ) . where the XDR ratios are high- 
est. Most of the absorbed photons in a PDR will heat 
the dust. An AGN (creating an XDR) generates a UV 
continuum which contains approximately ten times more 
energy than the X-ray field, and also heats the dust effi- 
ciently. In shocks, on the contrary, the gas is compressed 
and heated to higher temperatures, while the dust is not 
affected (except for shock velocities and densities that are 
orders of magnitude higher). Assuming that shocks are 
not heating the dust and that all the far-infrared lumi- 
nosity is reprocessed radiation from the AGN, we obtain 
a maximum AGN contribution of 10 — 15 percent. So, 
a shock dominated ISM can yield a much larger line-to- 
continuum ratio than PDRs and XDRs and this is exactly 
what we see in NGC 6240. 

Geometry of NGC 6240: The bulk of the gas mass does 
not coin cide with th e two AGN nuclei or with star for- 
mation. lEngel et al.l (|2010l ) relate gas masses, as traced 
by the CO J = 2 — 1 emission, to the different locations 
in NGC 6240: 5% and 25% to the northern and southern 
AGN nuclei, and 70% to the CO peak in the overlap re- 
gion, while using apertures of 1". Our FTS beam (> 19") 
is larger than the galaxy, and traces the total CO emis- 
sion. We determined a FWHM of 450 ±40 km s _1 for the 
CO J = 13—12 emission with a gaussianxsinc profile fit 
(larger than the instrumental resolution of 245 km s _1 



4 



TABLE 1 

Derived fluxes from the SPIRE FTS observations 

SUPPLEMENTED WITH GROUND-BASED OBSERVATIONS 



Line 




Z-obs 




^line 


-^linc 






[Ghz] 


[M m ] 


[Jy km s — 1 ] 


[L©] 


CO J = 1 - 


_ o a 




2600.76 


322 ± 29 


4.3 x 10 5 > d 


CO J = 2 - 


- l a 




1300.40 


1490 ± 250 


4.0 x 10 6 


CO J = 3 - 


- 2 a 




866.963 


3210 ± 640 


1.2 x 10 7 


CO J = 4 - 


-3 


449.99 


650.252 


4630 ± 370 


2.8 x 10 7 


CO J = 5 - 


-4 


562.41 


520.231 


5640 ± 150 


3.8 x 10 7 


CO J = 6 - 


- 5 


674.83 


433.556 


5910 ± 82 


4.8 x 10 7 


CO J = 7 - 


-6 


787.25 


371.650 


6010 ± 60 


5.6 x 10 7 


CO J = 8 - 


-7 


899.68 


325.225 


5830 ± 89 


6.3 x 10 7 


CO J = 9 - 


- 8 


1012.1 


289.120 


4770 ± 82 


5.7 x 10 7 


CO J = 10 


- 9 


1124.5 


260.240 


4160 ± 67 


5.6 x 10 7 


CO J = 11 


- 10 


1236.9 


236.613 


3160 ± 75 


4.7 x 10 7 


CO J = 12 


- 11 


1349.1 


216.927 


2590 ± 60 


4.2 x 10 7 


CO J = 13 


- 12 


1461.2 


200.272 


2080 ± 60 


3.6 x 10 7 


13 CO J = 


3-5 b 


645.21 


453.498 


< 112 ± 67 


8.6 x 10 5 


[Ci] 3 Pi - 


3 Po 


480.27 


609.135 


1750 ± 220 


1.0 x 10 7 


[Ci] 3 P 2 - 


3 Pi 


789.95 


370.415 


3330 ± 52 


3.1 x 10 7 


[Nil] 3 Pi - 


3 Po 


1426.1 


205.300 


3220 ± 77 


5.5 x 10 7 


H 2 2 U - 


2()2 


734.22 


398.643 


522 ± 52 


4.8 x 10 6 


H 2 2 02 - 


111 


964.61 


303.456 


798 ± 75 


9.2 x 10 6 


H 2 3i2 - 


303 


1071.5 


273.193 


548 ± 60 


7.0 x 10 6 


H 2 In - 


Ooo 


1086.7 


269.272 


234 ± 52 


3.0 x 10 6 


H 2 3 2i - 


3-12 


1135.3 


257.795 


606 ± 60 


8.2 x 10 6 


H 2 2 20 - 


2ii 


1199.8 


243.974 


394 ± 60 


5.6 x 10 6 


UID line 




1472.7 


198.689 


277 ± 60 


4.9 x 10 6 


H 2 v = 1 - 


S(l) c 




2.12 




3.5 x 10 7 


H 2 v = 1 - 


5(0) c 




2.22 




1.1 x 10 7 


H 2 v = 2- 


1 S(l) c 




2.25 




4.9 x 10 6 


H 2 v = 2 - 


1 5(2) c 




2.15 




2.2 x 10 6 



" IPapadopoulos et al.l l20TlT) 

upper limit 

ITecza ct al. (2000), extinction corrected values, luminosities ad- 
justed to adopted distance 

''relative uncertainties for Lii ne are the same as for Su ne 

at this frequency range). This is similar to the FWHM 
of the CO J = 2 — 1 line at the emission peak between 
the two nuclei, suggesting that the CO J = 13 — 12 emis- 
sion is also located in this region. The projected dis- 
tance between these nuclei on the sky is ~ 750 pc, and 
the true distance between them is estimated at 1.4 kpc 
(jTecza et al J 120 00). These authors assume that the two 
nuclei are on a circular orbit, of which the position angle 
and inclination are th e same as that of th e CO-disk be- 
tween the two nuclei (jTacconi et al.l[l999f ). and that the 
velocity difference between the two nuclei is 150 km s _1 . 
The AGN X-ray luminosities are not enough to power 
the CO excitation of the gas that is residing between 
the two nuclei, since they are geometric ally diluted and 
absorb ed: The luminosities derived by Komoss a et al.l 
(|2003|) give an X-r ay flux at 250 pc fro m the AGN of 
< 1 erg cm -2 s _1 . iVignati et al.l (|1999f ) derive a much 
larger luminosity by an obscured AGN, with an absorb- 
ing screen of Ah ~ 10 24 cm~ 2 which also reduces the 
X-ray fluxes below < 1 erg cm~ 2 s _1 . This leaves room 
for an XDR component located near the AGN nuclei, 
but even the most optimistic estimate of the X-ray lu- 
minosity is not enough to explain the combined H2 and 
CO luminosity of L = 2.5 x 10 9 L Q = 9.6 x 10 42 erg s _1 . 
This would require a strong coupling of the X-rays to the 
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Fig. 3. — C-shock model (solid line) with density njj = 5 X 
10 4 cm -3 and shock velocity v s = 10 km s _1 overlaid on the 
observed CO line fluxes (diamonds) for NGC 6240. 

molecular gas. 

3.1. Shock modeling 

Given our argum ents discussed above and the fact that 
vario us papers (cf. Jvan der Werf et al.lll99 3: Tccz a et al.l 
2000) have argued that C type shocks are exciting the 
H2 lines in NGC 6240, we analyze the CO ladder with 
shock models. Two types of shocks are used: a mag- 
netic continuous (C-type) and a non-magnetic jump (J- 
type) shock model. In both types, ro-vibrational lev- 
els of H2 are excited through high temperature H2— H2, 
H— H2, and He-H2 collisions. As the shock develops, 
the temperature of the gas becomes such that the ex- 
cited vibrational states become po pulated. We use the 
iFlower fc Pineau des Foretsl (|2003| ) shock code to model 
the H 2 v = 1-0 S(l)/v = 2 - 1 5(1) ratio and the 
CO ladder. IKristensen et al.l (|2007j) used this model to 
calculate a grid, spanning hydrogen number densities be- 
tween 71h = 10 4 and 10 7 cm -3 , velocities between v = 10 
and 50 km s _1 , and transverse magnetic held densities 
1/2 

b x n H /xG, with b=l and 5. The magnetic held den- 
sity relation implies values between 100 and 3000 /iG for 
the adopted densities, which is within t he range of valu es 
observed for galactic molecular clouds (|Crutcherl ll999). 

Using a rotational diagram and assuming that the CO 
lines are optically thin, we find T rot = 66 K for the 
CO rotational transitions between CO J = 5 — 4 and 
J = 7 — 6. This rotational temperature increases to val- 
ues T^t ~ 150 K at the highest CO transitions. We note 
that for the highest transitions the rotational temper- 
ature is a lower limit to the kinetic temperature, since 
these transitions are slightly sub-thermally excited. Us- 
ing a chi-square fit and including the errors provided in 
Table[U the best-fitting C-type shock model has a hydro- 
gen number density of nu = 5 x 10 4 cm -3 and velocity 
of v = 10 km s _1 (see Fig. [3]). The density is well con- 
strained, and chi-square values increase by an order of 
magnitude when going to densities n — 10 4 or 10 5 cm~ 2 . 
The uncertainty in the shock velocity is a few km s _1 . 
Downstream in this particular model, the CO gas has 
been compressed by a factor of 7.2 and the post-shock 
density of the CO emitting gas is 3.6 x 10 5 cm" 3 . 

The?; = 1-0 5(1) and v = 2-1 5(1) H 2 lines have up- 
per level energies of E = 6500 and 12500 K, respectively, 
and are therefore only excited when the gas temperature 
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T > 1000 K. It turns out that there is a degeneracy be- 
tween pre-shock density and velocity for reproducing the 
observed v = 1 - S(l)/v = 2-1 S(l) ratio of - 7 
(jTecza et a l. 2000) . The lower the pre-shock density, the 
higher the velocity that is required: The solutions range 
between hydrogen number density rtjj = 5 x 10 5 cm~ 3 
and a shock velocity v s — 47 kins" 1 to a pre-shock 
density of njj = 10 7 cm~ 3 combined with a velocity of 
v s = 16 km s . 

4. DISCUSSION AND CONCLUSIONS 

Shock modeling: Combining the model results for the 
H2 and CO emission, we find that the H2 v = 1 — 5*(1) 
to CO J = 10 — 9 intensity ratio is of order ~ 100. In 
contrast, the observed H2 / CO line luminosity ratio is 
approximately ~ 0.5. The low density, low velocity shock 
model fitting the CO lines, has a low temperature and 
does not produce H2 emission. From this we conclude 
that only a very small fraction of the gas mass is currently 
exposed to very powerful shocks (with either very high 
densities, jih = 10 7 cm -3 and moderate shock velocities 
v s = 16 kins" 1 or moderate densities, «h = 4xl0 4 cm' 3 , 
combined with a high shock velocity v s ~ 50 km s _1 ). 
Most of the shocked gas is settling and equilibrating with 
the ambient ISM, which is in agreement with the fast 
dissipation timescale derived below. 

Dissipation timescales: If we assume that all the gas 
is colliding at a shock velocity v s = 50 km s _1 (which is 
the highest velocity allowed by the models reproducing 
the H2 lines), the total amount of energy available (for 
an adopted gas mass of M = 1.6 x 10 10 M Q , in the 
middle of the LVG masses derived by various authors) is 
E = 0.5 M v 2 s — 4.0 x 10 56 erg, which is very similar to 
the value derived bv lTacconi et al.l (Qj)99). Assuming no 
additional energy input and that CO traces the bulk of 
the gas mass, this would imply that the CO shock energy 
is dissipated away within 6.6 million years, well within 
the orbital timescale of the two nuclei of 30 million years 
(|Tacconi et al.lll999D . 

Excitation by UV/X-rays vs. Shocks: The observed H2 
v = 1 - S(l) I v = 2 - 1 S(l) ratio falls within the 
range of ratios t hat are produced by X-ray dominated 
regions (see, e.g.. lMalonev et al"1ll996l Fig. 6). Also, the 
CO ladde r of NGC 6240 resembles the one observed for 
Mrk 231 (|van der Werf et al.l [20T0I) . However, as men- 
tioned before, NGC 6240 is lacking the bright OH + and 
H2O" 1 " lines, which are associated with gas clouds that 
are exposed to extremely high cos mic ray or X-ray ion- 
ization rates ([Mciicrink et al. 201]]), and the available X- 
ray photons are not sufficient to dominate the chemistry 
and thermal balance of the bulk of the gas (see Sect. 13. ip . 



The H2O lines are less luminous than in Mrk 231. This 
implies either lower water abundances or a less efficient 
mode of excitation. A full analysis is beyond the scope of 
this paper. These results and the similarity to the high- J 
CO line distribution in Mrk 231 suggest th at shocks pos- 
sibly due to the massive m o lecular outflow (iFischer et all 
1201(1 iFeruglio et al.l IMTct ISturm et al.ll201lD may also 
contribute to the CO line emission in Mrk 231. Indeed, 
base d on millimet e r inte rferometric observations of Mrk 
231, iCicone et all (|2012l ) note that in the inner region 
{R < 0.3 kpc), the CO(2 - l)/CO(l - 0) ratio is slightly 
higher, indicative of a shock contribution. 

Observing strategies for ALMA and other sub- 
millimeter facilities: Although the observed CO ladders 
for Mrk 231 and NGC 6240 are practically indistinguish- 
able, we argue that shock excitation and not X-rays are 
responsible for the excitation of the CO ladder, based on 
our knowledge of the geometry of the NGC 6240. Such 
an analysis will not be possible in the study of the ISM 
in galaxies at high redshift. At those distances, we are 
unable to resolve the offset between the two radio nu- 
clei and the location of the CO emitting gas, and X-ray 
fluxes cannot be determined. Therefore, we have to rely 
on dust continuum and line emission at far-infrared and 
sub-millimeter wavelengths. Here we have argued that a 
high line-to- continuum ratio is a key diagnostic for the 
presence of shocks. Additionally, more detailed mod- 
eling of emission by molecules and ions, such as H2O, 
H2O" 1 " and OFI + , will help in making this distinction, 
and are thus also highly recommended to include in fu- 
ture ALMA programs. 
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